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Abstract 

The effects of gravitomagnetic force on plasma oscillations are investi- 
gated using the kinetic theory of homogeneous electrically neutral plasma 
in the absence of external electric or magnetic field. The random phase as- 
sumption is employed neglecting the thermal motion of the electrons with 
respect to a fixed ion background. It is found that the gravitomagnetic 
force reduces the characteristic frequency of the plasma thus enhancing 
the refractive index of the medium. The estimates for the predicted effects 
are given for a typical white dwarf, pulsar, and neutron star. 

1 Introduction 

The outer crust of a compact star is formed of metal consisting almost exclu- 
sively of iron and other heavy elements. For a typical white dwarf the solid 
^^Fe core is surrounded by Si, O, C, and He ions and a gas of free electrons, 
whereas in neutron stars (pulsars), central densities are extremely high (up to 
10^^ kg /cm^) therefore the equation of state is relatively fixed. The outermost 
region of a neutron star is formed of a solid ^^Fe (density « 10kg /cm^) crust 
forming a crystalline lattice, surrounded by degenerate electrons. The star is 
thus regarded as enveloped in a highly ionized degenerate solid state plasma 
consisting of densely packed ions and relatively free electrons Tl. Collectively a 
plasma behaves as an oscillating system having a characteristic (resonant) fre- 
quency called the plasma frequency which for a degenerate solid state plasma 
is well defined throughout the medium. In particular the refractive index of a 
plasma medium is determined by the plasma frequency, which in turn deter- 
mines the transmission and reflection of radiation for the plasma On the 
other hand gravitational effects play a key role in all processes occurring in a 
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compact star, so that general relativistic effects are important for an adequate 
description of the phenomenon. Usually these effects are not directly observ- 
able but are manifest in an indirect way; for example via interaction with the 
magnetic field[3], in the accretion of matter^], and other material and radiative 
processes occurring in vicinity of the starjS). Since in these situations gravita- 
tional field of the star is comparatively weak (as compared to totally collapsed 
objects), and magnitude of the angular velocity of the star lies well below the rel- 
ativistic limit, the equations of the general theory of relativity can be expressed 
in an approximate form, adequate for a description of the gravitational field 
of the star. This approximation due to its formal analogy with classical elec- 
tromagnetic theory is called the gravitoelectromagnetic (GEM) approximation. 
The linearized geodesic equation contains a velocity dependent gravitomagnetic 
(GM) force whose origin lies in general relativistic frame dragging effects. It has 
been shown |H] that due to the gravitomagnetic force a coupling between the 
frequency of a material oscillator, consisting of a point mass, and the angular 
frequency of a gravitational source establishes. The coupling is especially en- 
hanced at the resonant frequency of the oscillator. It is therefore of theoretical 
as well as observational interest to study the possible effects of gravitomagnetic 
force on collective oscillations of a solid state plasma in compact star crustal 
regions where gravitational force is particularly effective. 

The paper is organized as follows. In the next section we explain the GEM 
approximation to the general theory of relativity and develop a force law, anal- 
ogous to classical electromagnetic theory, based on the geodesic equation. We 
then extend the formalism for a point mass oscillator, to include the collective 
behavior of the plasma in a GEM force field, using the kinetic theory for the 
plasma oscillations in section 3. The extension based on the random phase as- 
sumption leads to an analogous shift in the plasma frequency. Later in section 4 
we investigate the effect of the frequency shift on refractive index of the plasma 
as a single radially oscillating fluid. In section 5 we summarize the main results 
and discuss the possibility of detection of the predicted shift. Throughout we 
employ the gravitational units where G = 1 = c unless mentioned otherwise. 

2 The Gravitelectromagnetic Approximation 

In General Relativity the gravitational field of a massive object is described by 
the metric tensor gap- For sufficiently weak gravitating systems, such as the 
compact stars, a linearization the metric can be adequately assumed, where 

= Tlafj + hafS (1) 

and where rj^p = diag{-'\, —1, —1, 1) is the Minkowski metric tensor and hap is 
the perturbation to the metric such that hap <C 1 and = (x*, a;°) = (r,i) are 
the position coordinates of test particle, with time t being the affine parameter. 
Further requiring motion to be slow and F^^ denoting the Christoffel symbol, 
the geodesic equation 
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has the foUowing analogue in the GEM approximation[7l|S| 




(2) 



G + v X H 



(3) 



where 



G^-V(p, H = Vx4a 



(4) 



and 



if = 





(5) 



p being mass density and V is the volume. Remarkably the formal analogy of the 
above results with the classical electromagnetism extends to the a gravitational 
field equations. 

Since in the slow rotation approximation the deformations in mass distri- 
bution in the star due to rotation are negligible, the star can be regarded as 
a slowly rotating sphere of homogeneous mass density M . Then the gravito- 
electric (GE) force acting on a unit mass in vicinity of the star is given by the 
Newtonian gravitational force 



and H is given by 

ll = -fMR^in.rI,-l^) (7) 

where R is the radius and is the angular velocity of the gravitational source. 

The second part of equation (3) has a non-Newtonian origin and is regarded 
to exhibit typically general relativistic effects, such as the Lense-Thirring drag- 
ging of frames. Within the GEM approximation the effect can be interpreted as 
'gravitomagnetic current' induced in the vicinity of the gravitational source due 
to its rotation. It therefore plays an important role in testing Einstein's theory 
of gravitation in the weak field and slow rotation approximation. Further it 
has been demonstrated^ that the form of the GM potential does not depend 
on the choice of a particular frame or coordinate system used, thus making it 
physically significant for all material as well as radiative processes. However to 
measure the effects of this force on a given physical system high accuracy in 
experiments is required 10 . 
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3 Equation of Motion for Collective Plasma Os- 
cillations 



A detailed kinetic theory of plasma oscillations in metals has been developed 
using the Fermi-Dirac statistics leading to the quantized theory of elementary 
excitations 11 . However the main features of collective plasma oscillations due 
to the effects of the gravitomagnetic field can be discussed within the classical 
treatment. We consider the phenomenon of radial electron oscillations confined 
to two dimensions, whereas the angular velocity vector of the star is inclined 
at an angle x to the plane of oscillations. The components of acceleration due to 
gravitoelectric force G to be of constant magnitude g and the gravitomagnetic 
force on the ith electron has the magnitude fiHsinx- The equations of motion 
for the radial oscillations of the ith electron can thus be written as 



y^{\)exp[ik{ri - r^)] - g + ViHsinx- (8) 
m ^-^ k 

j-k' 

where e is the electronic charge and m the electron mass. The prime over k 
denotes a sum where fc = is excluded. Since acceleration due to gravitoelectric 
force is of a constant magnitude, its effect on the electron oscillations with 
respect to the fixed background of ions will be negligible. We therefore redefine 
fj as fj — g in the equation of motion (8). For all practical purposes electrons 
can be considered as point particles, therefore in an arbitrary region of unit 
volume the particle density is given bv|12| 

p(r) = ^<5(r-rO, (9) 

i 

We decompose the particle density function into Fourier components as 

/pWe.p|-.t.l* = i:e.p|-..nl. (10) 

i 

such that 

p{r) — exp[ik{r — r^)] (11) 

Here it is understood that i as a multiplier is the complex number denoting 
and not the script employed for the summation. Differentiating Pi^{ri) 
with respect to time we obtain 

Pkin) = -i^{kf,)exp[-ikn], (12) 

i 

and 

Pki'T'i) = - ^[(kh)'^ + ikri]exp[-ikri]. (13) 
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A substitution from (8) into (13) gives 



Pkin) = - ^{{kr-if + ik{ — '^{—)pkexp[ik'ri] + fiHsinx)}exp[-ikri 



(14) 

which can be simpUfied to give 



Pkin) = - y](fcri)^ea;p[-ifcri] - ^{^){exp[i{k' - k)ri]}exp[-ik'rj] 

% ijk' 

—i kf-jH sinxexp[—ikri] . (15) 

i 

Here the second term on the right hand side can be split into two parts one for 
k' = k and the other for k' ^k: 

^^^^ exp[—ikri] + A exp[—ik'ri\ exp[i{k' — k)rj]. (16) 

The phase factors, randomly distributed in the complex plane, approximately 
cancel each other out. Therefore the sum of such factors is small, and hence the 
second term being a product of such summations is negligible. Thus random 
phase approximation gives in equation (15) 

Pk{ri) = —'y^^{kfi)^exp[—ikri] — ^^^'^ ea:p[— ifcrj] 

i i 

—i kriexp[—ikri]Hsinx- (17) 

i 

Further since in a metal electrons are strongly bound to the ions, their 
random thermal motion is negligible, therefore the contributions coming from 
the first term containing fc^ arc also negligible. Using equations (10) and (12) 
in the second and third term of equation (17) we obtain 

hin) = —Pk + hHsinx- (18) 



4 Plasma Frequency Shift and Refractive Index 

To determine the resonant frequency for the Fourier component for electron 
density function Pki^i) we assume real functional dependence on time t in the 
form of the solution Akcosuit to equation (18), where is the amplitude of the 
density fluctuation and u is applied external frequency. The substitution leads 
to the following condition for resonance 
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ojI-lu^ - ojHsinx = 0. (19) 

Since H <C ujp, solving (19) with uj unknown and neglecting terms involving 
squares and higher powers of H/wp, we obtain for the plasma frequency 

H , , 

w = Wp — ^sinx- (20) 

Clearly the shift is dependent upon the rotational frequency and mass of the 
star, and also upon the angle of inclination x- It follows from expression (7) 
that at the surface of the star r = R, the gravitomagnetic force has magnitude 
given by 



H=\U\c^ /iVl + 3cos2 X, (21) 

where fj, = {AGM/^R(?)Q., Q. being the magnitude of the angular velocity vector. 
Substituting from expression (21) into (20) we obtain an expression relating 
plasma frequency ujp to the angle 

— ^l-^^Vl + 3cos2x, li^ujp. (22) 

U3p Zujp 

Here the requirements H <^ cOp and fi ^ ojp arc generally valid for typical 
compact stellar sources with ^ ranging from 0.1687 x 10~^ Hz (for a typical 
white dwarf of mass IMg, = 1.989 x 10^^ kg , radius 7 x lO^rn and angular 
frequency li? 2;) to 236.29327Jz (for a typical neutron star of mass 2Mg,, radius 
1 x lO^m, and angular frequency IkHz) with corresponding plasma frequency 
ranging approximately 5.65 x lO^Hz to 5.65 x lO^Hz or above. Denoting cOp = 
ujp — H/2, where H is given by (21), we have in terms of plasma frequency we 
have the following expression for the index of refraction e of the medium 



.(X) = Ji - i^r = J^-i^- f^^ATw^)^ (23) 

where LUa is the angular frequency of waves (e.g. electromagnetic) propagat- 
ing through the plasma; whose transmission or reflection, for a given compact 
source, is determined by the usual conditions 

i^a < i^p, reflection; (24) 



LJa > i^p, transmission. (25) 

The shift in plasma frequency in the equatorial plane (x = 7r/2) is ujp — ii/2 
whereas in the plane orthogonal to it (x = 0) the shift is Wp — ji. Denoting 
the frequency shift in the equatorial plane by Cjp\\ and by u}p±_ for the orthogo- 
nal plane, we find that the maximum difference in the frequencies transmitted 
through the plasma: 

I ^pl. - ^p\\ |max= ^- (26) 
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In figures (1) and (2) we give plots for the typical cases of neutron star, 
pulsar, and white dwarf, between the plasma frequency shift Lo/ujp and the 
angle x (using expression (22)) and between the plasma refractive index e as a 
function of the angle x (using expression (23)) respectively. 

5 Conclusions 

In view of the above investigations we find that the 'gravitomagnetic oscillations' 
have a coupling to the plasma oscillations which become especially effective at 
the plasma frequency. This is so even if the plasma co-rotates with the star 
mainly because of the fact, that the gravitomagnetic force as well as the elec- 
tron oscillations are along the radial direction, as is evident from expression 
(3) and (8). It is clear from figure (1) and (2) that the plasma frequency the 
gravitomagnetic force produces a reduction in the frequency, thus enhancing 
the refractive index of the plasma. Further the dependence of the shift on com- 
ponent of the angular frequency vector of the gravitational source to the plane 
of oscillations of the plasma gives a difference in the refractive index, which 
results in changing the allowed range of frequencies transmitted through the 
plasma, by a maximum amount /x/2. For a typical white dwarf /i is approxi- 
mately 0.1687 X lQ~^Hz^ for a typical pulsar \.QhAOHz and for a neutron star 
236.2932ifz. The predicted difference in the frequency of radiation emitted by 
a compact star though still very small, may possibly be observed, for example in 
radiation spectra of various compact sources, as the observations become more 
refined. 
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Figure Captions: 

Figurel: Plots for the shift oj/iOp in the plasma frequency of a compact 
star atmosphere as a function of the angle of inclination x of the plane of 
observation to the angular frequency vector for the case of a white dwarf 
(M = IM® = 1.989 X lO^o/cg, i? = 7 x lO^m, VL = IHz, Wp ^.65 x lO^Hz, 
a pulsar {M = 1.4M®, i? = 3 x lO^'m, fl = 30Hz, Up = 5.65 x lO'^ H z) ,a.iid 
a neutron star (M = 2M@, R=lx lO^'m, n = IkHz, ujp = 5.65 x l(fHz) 

Figure 2: Plots for the refractive index e of a compact star atmosphere 
as a function of tlic^ angle of inclination x of th'- plane of observation to 
the angular frequency vector for the case of a white dwarf (M = IMg, = 
1.989 X 10^°%, i? = 7 X lO^m, Q. = IHz, tOp =6.65 x lO^i?^, a pulsar 
(M = 1.4M®, i? = 3 X 10'*m, SI = 30Hz, ujp = 5.65 x 10'^Hz),and a 
neutron star (M = 2M@, i? = 1 x lO^m, n = IkHz, iOp = 5.65 x lO^ifz) for 
the propagation of an electromagnetic signal of angular frequency IQ^Hz, 
IQ^Hz, and lO^Hz respectively. 
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